NMR as a probe of the relaxation of the magnetization in magnetic molecules 
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We investigate the time autocorrelation ol the molecular magnetization M(t) for three classes of 
magnetic molecules (antiferromagnetic rings, grids and nanomagnets) , in contact with the phonon 
heat bath. For all three classes, we find that the exponential decay of the fluctuations of M(t), 
associated with the irreversible exchange of energy with the heat bath, is characterized by a single 
characteristic time t(T, B) for not too high temperature T and field B. This is reflected in a nearly 
single-lorentzian shape of the spectral density of the fluctuations. 

We show that such fluctuations are effectively probed by NMR, and that our theory explains the 
recent phenomenological observation by Baek et al. (PRB70, 134434) that the Larmor-frequency 
dependence of 1/Ti data in a large number of AFM rings fits to a single-lorentzian form. 
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PACS numbers: 76.60.Es,76.60.-k,75.60. Jk 



A central aspect of the physics of magnetic molecules 
is the way molecular observables, and particularly the 
magnetization M, are affected by interactions of the spins 
with other degrees of freedom. When the latter behave as 
a heat bath, they cause decoherence of the time evolution 
of molecular observables, usually leading to relaxation 
dynamics and exponential time-decay of equilibrium fluc- 
tuations. The understanding and characterization of 
relaxation mechanisms in crystals containing magnetic 
molecules is not only important for evident fundamental 
reasons, but also because relaxation and, more generally, 
decoherence, constitute a major obstacle in the envisaged 
technological applications of these molecules. In partic- 
ular, relaxation of M through phonons sets an upper 
bound lj to the temperature range in which nanomagnets 
could be used as classical bits in a memory. In the pos- 
sible application of some molecules as elementary units 
(qubits) in quantum information processing (QIP)^, en- 
ergy dissipation from the qubits to the nuclear-spins heat 
bath is the main concern regarding the efficiency of quan- 
tum logical gates and the performance of the computa- 
tion. 

In this work, we show that for three important classes 
of magnetic molecules (antiferromagnetic (AFM) rings, 
grids and nanomagnets), the exponential decay of the 
fluctuations of M(t) due to spin-phonon interactions is 
dominated by a single characteristic time t(T, B) over 
a wide range of temperatures and applied fields. While 
this was already known for the long-time decay of M 
in nanomagnets, our results show that this applies to 
completely different classes of molecules, and holds on 
much shorter time-scales as well. Indeed, in the consid- 
ered (T, B) regimes this single characteristic frequency 
Xq = 1/t dominates the fluctuations spectrum to such a 
large extent that even experimental probes whose intrin- 
sic frequency is very different from Ao will nevertheless 
detect mostly this single spectral contribution. In partic- 
ular, we show that NMR measurements of the longitudi- 
nal nuclear relaxation rate 1/Ti can be used as a direct 
probe of the fluctuations of M, and that experimental 



1/Ti data are fully consistent with our picture. For in- 
stance, we explain the recent experimental observation 
3] that the Larmor-frequency dependence of 1/Ti data 
in a large number of AFM rings fits to a lorentzian form, 
and we elucidate the meaning and behavior of the char- 
acteristic time r(T, B) of this lorentzian. In addition, we 
collected new data on the Mn3x3 grid which give further 
support to our picture. At last, we reanalyzed NMR data 
on the prototype nanomagnet Fes- 

Each magnetic molecule can be described by the spin 
Hamiltonian 



i>j=l,N i=l,N k,q 

^ s, ■ Dy ■ s 3 - gfi B BS z , 

i>j=l,JV 

Ih 



(1) 



where Sj are spin operators of the i magnetic ion in 
the molecule. The first term is the isotropic Heisenberg 
exchange interaction. The second term describes local 
crystal-fields (CFs), with Oj?(si) Stevens operator equiv- 
alents for the z-th ion£| and bt(i) CF parameters. The 
third term represents dipolar anisotropic intra-cluster 
spin-spin interactions. The last term is the Zeeman cou- 
pling with an external field B with S z = X)i s i,z- The 
Heisenberg contribution is usually largely dominant, and 
therefore the energy spectrum of H consists of a series of 
level multiplets with an almost definite value of the total 
spin S. 

Much of the physics of nanomagnets such as Fes can be 
described to a good degree of approximation by retaining 
the lowest 5-multiplet (e.g., S = 10 in Fe 8 )0. In AFM 
rings such as Fe6 and Cr 8 and in grids, several distinct S- 
multiplets must be included for a good description. The 
ground and first-excited multiplets have S — and 5=1 
in AFM rings 0, while they have 5" = 5/2 and 5" = 7/2 
in the Mn3x3 gr id0. In rings and grids the anisotropic 
terms in Eq. (1) have no qualitative effects for most ob- 
servables, (including those on which we focus here) and 
in most situations (i.e. away from intermultiplet (anti)- 
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crossings caused by B). In view of this, and since numer- 
ical calculations based on the full Hamiltonian Eq. (1) 
are cumbersome, in these two classes of systems we will 
use Eq. (1) neglecting the anisotropic terms. In addi- 
tion, large-scale numerical techniques must be used for 
Mn3x3, resulting in a limitation in the range of T which 
can be studied. 

Our first goal is the calculation of the damping of equilib- 
rium fluctuations of the molecule magnetization M(t) = 
gS z due to interactions of the electron spins with a heat 
bath. We adopt the well-established framework pre- 
sented in Q for the irreversible evolution of the density 
matrix p s t(t) (in the representation of the exact eigen- 
states of (1)) in presence of a heat bath. Here we focus 
on time scales t much longer than the typical periods 
of free evolution of the system (2wh/(E s — E t ), with E s 
and E t eigenvalues of (1)). In this case, the so-called 
"secular approximation" enables the time evolution of 
the off-diagonal matrix elements p s /t (t) to be decoupled 
from that of the diagonal elements p ss (t) :— p s (t). In a 
frequency-domain picture (e.g., the outcome of a scatter- 
ing experiment), this corresponds to a clear separation 
of inelastic (IE) and quasielastic (QE) spectral contribu- 
tions. Then the (single-lorentzian) broadening of IE tran- 
sitions merely reflects the lifetimes of the involved pair 
of levels (basically determined by Fermi's golden rule ap- 
plied to the spin-bath interaction). These IE off-diagonal 
contributions to the evolution equations for M(t) will be 
neglected since we are here interested in the subtler QE 
components of the dynamics. Only the latter are de- 
tected by low-frequency probes such as AC-susceptibility 
or NMR. 

In the secular approximation the populations of the 
molecular eigenstates, p s (t), evolve through the well- 
known rate (master) equations 

p s (t) = ^2W stPt (t), (2) 
t 

where W s t is the st element of the rate matrix W, i.e. 
the probability per unit time that a transition between 
levels \t) and |s) is induced by the interaction with the 
heat bath (W ss =-J2t W ts)- 

We are here focusing on temperatures for which relax- 
ation is determined by interactions of the spin degrees 
of freedom of each molecule with the crystal phonon 
heat bath. Since we are interested in general quali- 
tative features of the relaxation dynamics and not in 
compound-specific details, we will eventually assume a 
unique simple form for the spin-phonon coupling for all 
compounds. The main contribution to this coupling 
comes from the modulation of the local rank-2 crystal 
fields ^(O^f (si) by phonons : 

V= £ E^Qr^Cr^, (3) 

i=Z,Arri, 7i 

where Qr, , 7i are symmetrized displacements of neighbor- 
ing atoms of the Zth magnetic ion, labelled by irreducible 



representations Ti of the local point group of ion I, 7 ; 
being the component. Cr i!7i are linear combinations of 
0|(s;) transforming as 1^,7;. For instance, in octahe- 
dral local environment, if T = r$, Cr,i = Oj^s;)/-^ and 
Cr,2 = 02 (s;). The symmetrized displacements Qr^-yi 
of the neighboring atoms can be expanded in terms of 
normal coordinates qk a of the full crystal lattice, 

Qr,,n = ^a^a{l)q-ka, (4) 

ko- 

where a-^ a (l) are known as Van Vleck coefficients. This 
leads to the following general form for the spin-phonon 
coupling potential 

v = E E E^'k-^^^+cUj, (5) 

l=l,NQ=-2,2 kcr 

where Cq(1, k, a) is the coupling constant between the Q- 
type electric quadrupole on ion I (0^(s/)), and phonon 
modes of wavevector k and branch index ct|9(. Clearly, 
a detailed microscopic calculation of Cg(Z,k, a) is un- 
feasible: it would require to model the (complicated) 
phonon spectrum, and to fix the value of a great number 
of parameters Vy i . Yet, experimental information is com- 
pletely insufficient to accomplish this. Therefore, we have 
chosen as usual^tj to adopt a Debye model for phonons, 
and for the sake of simplicity we assume Q = 0,1,2 with 
0Q(Z,k) independent of I and Q, although compound- 
specific Cq(1, k, a) would be needed for precise quantita- 
tive calculations, which are beyond the aim of this work. 
We have checked that the results do not change qualita- 
tively if different choices are made |l l|. 
The rate matrix W appearing in Eq. (2) can be cal- 
culated by first-order perturbation theory applied to V, 
leading to 

^ st = 7 7r|(i|E E 2 Q (s,)k)| 2 A^(A st ), (6) 

i=l,N Q=0,l,2 

with n(» = (e^ Hx - 1) -1 , A st = (E s -E t )/h. In our sim- 
plified spin-phonon model 7 is the unique free parameter 
and describes the spin-phonon coupling strength. It will 
be determined by comparison with experimental data. 
The next step is the calculation of the Fourier-transform 
Sa,b(w) of the cross-correlation function (AA(t)AB(Q)) 
implied by Eq. (2) for equilibrium fluctuations of two 
generic observables A, B. We exploit the fluctuation- 
dissipation theorem, which implies 

Sa,b M = = 2k B TR' A . B M , (7) 

with the imaginary part of the susceptibility, 

and R'ab{z) the real part of the Laplace transform of 
the relaxation function R.A,B(t)- The latter is defined 
as \im e _ ((A(t)) — (A) eg )/e, with a perturbation eB 
abruptly switched-off at t = 0. We have calculated 



3 




5 1 1 5 2 2 5 3 ~ 0.5 1 1.5 2 2.5 3 




05 0.1 0.15 2 25 3 35 1 2 3 



FIG. 1: Calculated frequency weights A(Xi, T, B) of the mag- 
netization autocorrelation vs A/T (a,b), with A the exchange 
gap, or vs 1/T (c) (x-axis). The j/-axis is log 10 (A) (in THz). 
The color maps log 10 (4(A 1 ,T,B)/(xT)L(see Eq. (9)). The 
hamiltonian parameters are given in (alHQ- F° r each value 
of T, the spectra have been normalized by \T, which is pro- 
portional to the size of equilibrium fluctuations. The white 
dashed lines are the 7 Li (a), 35 C1 (b), and X H (c) L armor an- 
gular frequencies. The results for Cr8 (not shown) are similar 
to those of Fe 6 Li, with 7 = 5 x l(T r Thz~ 2 . 
Panel (d) shows the characteristic frequency Ao vs A/T for 
Cr8, FeeLi, FeeNa and Mn3x3 extracted from the experimen- 
tal 1/Ti data (see text). The data are taken from ref. Q, 
apart from those for Mn3x3 which are new. The red line 
represents the dominant frequency Ao of FeeLi deduced from 
panel (a). 



RA,e(t) by solving the system of equations (2), obtaining 
the following expression: 

q,t 

x(A qq - (A)M)Re (pi e "(- 1 ¥ )^) 1 (8) 

where p[ is the equilibrium population of the t-th level, 
andStt = (t\B\t). This leads to correlations decaying as a 
sum of exponentials with n characteristic times Tq E given 
by the inverse eigenvalues of — W, n being the dimen- 
sion of the molecule spin Hilbert space. Therefore, the 
QE frequency-spectrum is in general a sum of lorentzians 
centered at zero frequency having width 1 / Tq E . It should 

(i) 

be stressed that the n Tq E cannot be directly determined 
from the n level lifetimes T^h. Although both Tq E and 

(s) 

T iife con tain integrated information on the spin-bath dis- 
sipation channels, this information is different for the 
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FIG. 2: Up : calculated 1/Tj./(xT) as a function of the 
temperature T for Org, Fee Li and Mn3x3 for different fields 
and nuclei. Inset : Calculated l/Ti(T) (in a.u.) for Feg 
for different values of B applied along the easy axis. Down: 
Measured l/T x /(xT) vs T for Cr 8 , Fe 6 Li and Mn3x3. In order 
to fit the figure, data for FegLi were divided by 3, and those 
for Mn3x3 by 5. 

two sets. Thus, low-frequency measurements probing the 
QE dynamics provide different and complementary infor- 
mation with respect to direct lifetime measurements by 
high-frequency techniques (e.g., inelastic neutron scatter- 
ing). 

We focus now on the important case A = B = S Zl i.e. 
we study the spectrum of fluctuations of the cluster mag- 
netization. In general, many different may be ex- 
pected to contribute to this spectrum, i.e. 

S a „s,(u,T,B)= M*i,T,B) x( ^: 2 B l 2 , (9) 

%— l,n v 7 

where A; (T, B) = 1 /tq B . Figure 1 shows intensity plots 
of the frequency weights A(Xi, T, B) of the magnetization 
autocorrelation, calculated from Eq. (8), as function of 
temperature T for three representative systems. The sin- 
gle free parameter 7 has been estimated by the observed 
peak in the T-dependence of the NMR 1 /T\ (see below) , 
apart from Fes f° r which relaxation data allow to fix 7 
directly. The precise value of 7 does not change the spec- 
tra appreciably, but merely sets the frequency scale. 

The calculations reported in Fig. 1 show that for a 
fairly wide range of values of T and B, a single frequency 
Xq(T,B) (and therefore a single lorentzian) dominates 
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the spectrum in the three classes of molecules we con- 
sidered (AFM rings, AFM grids and nanomagnets). The 
dominating frequency Ao is different from the lifetime 
broadening of any level. In the special case of nanomag- 
nets (Fes, M1112), the calculated Ao coincides with the 
smallest A^. Hence, in this case Aq 1 is also the time de- 
termining the asymptotic return to equilibrium of the 
magnetization when the nanomagnet is prepared in a 
highly nonequilibrium state by a strong applied field. In- 
stead, our calculations show that in rings and grids Ao 
is orders of magnitude larger than the smallest Aj (the 
latter is not visible in Fig. 1 since the associated weight 
is tiny). AFM rings are characterized by a single domi- 
nating characteristic frequency Ao for B <C B c with B c 
the field at which an S = 1 level crosses S — 0, and for 
T < 12J/N. This corresponds to the regime in which the 
magnetization and its fluctuations are mainly determined 
by levels belonging to the S = 1 manifold. In this single- 
frequency regime, the decay of fluctuations takes places 
mostly through indirect Orbach-like processes involving 
higher-lyings S'-manifolds, resulting in an approximately 
exponential dependence in 1/T of Ao (see Fig. 1). A dif- 
ferent T-dependence occurs at low T, where direct tran- 
sitions within the field-split S = 1 manifold constitute 
the main relaxation mechanism. Fes is also character- 
ized by a single dominating characteristic frequency Ao, 
depending exponentially on 1/T, if the Zeeman energy is 
small with respect to the anisotropy barrier. Finally, the 
Mn3x3 grid is again characterized by a dominating fre- 
quency for B <C B c with B c the field at which an S = 7/2 
level crosses S = 5/2, and for T < A/k B , with A - 0.9 
meV the exchange gap. Therefore, for all these systems 
the calculated dynamical structure factor of Eq. (9) is 
dominated by a single contribution: 

S St ,s. (a;, T, B) ~ A(A 0) T, B) - ( *° ^ ^ . (10) 

From linear-response theory A(\q,T, B) oc x(B,T)T, 
with x(B,T) the static Curie susceptibility. Therefore, a 
measurement of Ss z ,s z allows the T- and B-dependence 
of Ao(T, B) to be extracted straightforwardly. In the fol- 
lowing we will show that NMR confirms this picture and 
allows Ao(T, B) to be determined. 

The nuclear spin-lattice relaxation rate 1/Ti can be cal- 
culated by a perturbative treatment of the hyperfine 
interactions between nuclear and electronic spins [l2^. 
Within this framework, 1/Ti is given by a linear com- 
bination of the Fourier transforms of electronic spin cor- 
relation functions S Si atS . (3 , evaluated at the Larmor fre- 
quency lul(B). By considering only the Heisenberg and 
Zeeman terms in H (Eq. (1)), the spin lattice relaxation 
rate of a magnetic nucleus of Larmor frequency lol inter- 
acting with a cluster composed of N spins may be written 
as 

1/Ti cx ^j(SsI^(ujl) 

i,j = l,N 



+S s *,s*(-vl)) + Pij{S s + .-(wi) + S + -(-ul) 

1 3 i ' 3 * ' 3 

+S s - s +(lo l ) + S s7s +(-u l )), (11) 

where oiy and /3y are geometrical coefficients of the 
dipolar interaction between nuclear and electronic spins. 
Since generally u>l <§; A st , only the QE part of the 
Sa,b(w) in the above formula contributes to 1/T%. If 
anisotropy is neglected, Eq. (8) shows that only the 
a = 13 = z terms are nonzero because (t\si^ x \t) — 
(t\si tV \t) = 0. In addition, for a magnetic ring (£|sf |£) = 
(t\sj\t) = l/N(t\S z \i), where S z is the z component of 
the total spin operator S. Then, 

l/T 1 = GS Sz ,s z (u L ), (12) 

where G is a scale factor depending only on the positions 
and spins of the probed nuclei. The same result holds 
also for grids as long as the ground multiplet gives the 
main contribution (i.e. at sufficiently low T). In fact, in 
this regime all S Si z<s . t (u>) oc Ss z ,s z ( w )- We have numeri- 
cally checked that this holds approximately in the single- 
lorentzian regime (see Fig. 1). If magnetic anisotropy 
cannot be neglected, as in the case of Fes, the simple 
result of Eq. (12) can be recovered if B is along the 
easy axis, by exploiting the fact that within the ground 
S'-multiplet S Si Sj (lj) oc Ss z ,s z ( UJ )- Therefore, in all 
these systems NMR measurements of 1/Ti can be used 
as a probe of Ss z ,s z ( UJ ). 

Eqs.(10) and (12) imply that, as a function of T, 
l/Ti/(xT) should display a single peak when A(T, B) — 
u>l(B). Moreover, since X(T,B) depends little on B for 
not too large B, l/T\/{x7) measured on different nuclei 
and by setting the value of B in order to keep ojl(B) 
fixed, should display a peak at approximately the same 
T . This is clearly visible, for instance, in our calcula- 
tions of l/Ti/(xT) for FegLi on Li and H nuclei, shown 
in Fig. 2, which contains our calculated T-dependence of 
1 jT\ . By increasing B by a factor p and keeping the nu- 
cleus fixed, the peak value of l/T\/{yT) should decrease 
by about p. This can be seen in Fig. 2 in the calculated 
curves for Crs. This is what had been experimentally ob- 
served in , where a phenomenological single- lorentzian 
expression for the Larmor- frequency dependence of 1/T± 
had been shown to fit NMR data in several AFM rings 
very well. Part of the data are included in the lower 
panel of Fig. 2 together with new data we collected on 
the Mn3x3 grid. The agreement with our calculations 
is very good. In addition, panel (d) in Fig. 1 contains 
the T-dependence of the characteristic frequency Ao ex- 
tracted from the lorentzian fit of the data|3j. The results 
are in good agreement with the behavior of the calcu- 
lated Ao in panels (a,b). Similar agreement is obtained 
for the Crs ring. We remark that in the T-range where 
the peak in 1/Ti occurs the calculated T-dependence of 
Ao is exponential in 1/T (relaxation through intermul- 
tiplet Orbach processes), and in good agreement with 
experimental data. Hence, the phenomenological power- 
law T-dependence assumed in 0], although fitting the 
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data over a limited temperature range, is not confirmed 
by our microscopic calculation. 

In the case of Fes, the free parameter 7 has been fixed in 
order to reproduce the T-dependence of the relaxation 
time of the magnetization in the thermally-activated 
regime measured in There is therefore no free pa- 

rameter in the calculation. As can be seen in Fig. 1, the 
calculated Ao(T) is much smaller than the 1 H Larmor 
frequency. Indeed, measured 1/Ti data do not show a 
peakH, and are in good agreement with our parameter- 
free calculations reported in the inset of Fig. 2. 
In conclusion, a general expression for the quasielastic 
part of the dynamic structure factor in presence of a heat 
bath has been derived, starting from the master equation 
applied to the full microscopic hamiltonian. By special- 
izing this expression to the case of the total molecule 
magnetization, we have calculated (Fig. 1) the frequency 
spectral weights of the decay of fluctuations of this ob- 



servable. We have shown that in AFM rings, grids and 
nanomagnets the decay of equilibrium fluctuations due 
to interactions with phonons is characterized by a sin- 
gle characteristic frequency Xq(T,B) over a wide range 
of temperatures and fields. We have also shown that the 
NMR spin-lattice relaxation rate 1/Ti can be used to 
measure Xo(T,B). This microscopic result differs from 
the previous phenomenological models of 1 /T\ in rings 3J 
and nanomag nets [3, where 1/Ti was assumed to probe 
level lifetimes in presence of phonons, and not a relax- 
ation time. Our calculations are in good agreement with 
existing experimental NMR data, and with our new mea- 
surements on the Mn3x3 grid. 

Useful discussions with F. Borsa are gratefully acknowl- 
edged. We thank L.K. Thompson for the synthesis of 
Mn3x3, for which preliminary NMR results are published 
here. 
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